Symmetries and the conductance of graphene nanoribbons with long-range disorder 
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We study the conductance of graphene nanoribbons with long-range disorder. Due to the absence 
of intervalley scattering from the disorder potential, time-reversal symmetry (TRS) can be effectively 
broken even without a magnetic field, depending on the type of ribbon edge. Even though armchair 
edges generally mix valleys, we show that metallic armchair nanoribbons possess a hidden pseu- 
dovalley structure and effectively broken TRS. In contrast, semiconducting armchair nanoribbons 
inevitably mix valleys and restore TRS. As a result, in strong disorder metallic armchair ribbons 
exhibit a perfectly conducting channel, but semiconducting armchair ribbons ordinary localization. 
TRS is also effectively broken in zigzag nanoribbons in the absence of valley mixing. However, we 
show that intervalley scattering in zigzag ribbons is significantly enhanced and TRS is restored even 
for smooth disorder, if the Fermi energy is smaller than the potential amplitude. The symmetry 
properties of disordered nanoribbons are also reflected in their conductance in the diffusive regime. 
In particular, we find suppression of weak localization and an enhancement of conductance fluctu- 
ations in metallic armchair and zigzag ribbons without valley mixing. In contrast, semiconducting 
armchair and zigzag ribbons with valley mixing exhibit weak localization behavior. 

PACS numbers: 73.63.Nm, 72.80.Vp, 73.20.Fz, 73.23.-b 



I. INTRODUCTION 

The bulk electronic properties of graphene^ are modi- 
fied significantly in nanoscopic samples, where the influ- 
ence of the edges becomes decisive. Edge effects are most 
prominent in narrow graphene nanoribbons (GNRs) with 
the boundary structure determining whether the elec- 
tronic spectrum is semiconducting (gapped) or metallic 
(gapless) PHU The nanoribbons with the highest symme- 
try with regard to the graphene lattice exhibit zigzag 
and armchair edges, shown in Fig. [T] While zigzag 
GNRs are always metallic, armchair GNRs are catego- 
rized into 'metallic' or 'semiconducting' depending on 
their widthPUl GNRs of an orientation in between arm- 
chair and zigza g ha ve been shown to effectively behave 
as zigzag GNRs™ 

The first experiments on lithographically defined 
GNRs failed to find conclusive evidence for these edge 
effects,™ but since then great effort has been spent on 
improving the edges of GNRs: scanning tunneling mi- 
croscopy has been used to tailor edges™ GNRs have 
been chemically derived from solution phased they have 
been obtained by 'unzipping' of carbon nanotubesp^' 
cut out of graphene sheets by anisotropic etching using 
nickel clusters^ or by sonochemically methods,^ or they 
were self-assembled inside carbon nanotubesPO (For an 
extended overview we refer to a recent reviewP^) 

In certain situations, for example when the dynamics 
of a system is chaotic or diffusive, its quantum transport 
properties are mainly determined by very few symmetries 
of the system, in particular the presence or absence of 
time-reversal symmetry (TRS) TRS is usually broken 
by magnetic fields. However, Berry and MondragorP^J 
showed that in a (hypothetical) neutrino billiard, TRS is 
broken even in the absence of a magnetic field. In fact, in 



this case TRS is broken by the edge of the billiard itself. 
The Dirac equation used in Ref. [19] to model neutrinos 
corresponds exactly to the effective low-energy Hamilto- 
nian of a single valley of graphene, prompting efforts to 
realize such an effective TRS breaki ng in graphene, for 
example by different kinds of disorder^ 20 * 21 * or edgesJ^f 

Zigzag GNRs have also been shown to exhibit this kind 
of effective TRS breaking when only long-range disor- 
der is present, such that the valleys remain uncoupled 
and a single-valley Dirac equation description is appli- 
cable. In this case, (effective) TRS is broken by the 
zigzag edges, placing zigzag GNRs with long-range dis- 
order into the unitary symmetry class (no TRS) The 
symmetry class was also shown to influence the conduc- 
tance of zigzag GNRs dramatically: For long-range disor- 
der zigzag GNRs exhibit a perfectly conducting channel 
(PCC), i.e. a minimum of one conductance quantum even 
in the strongly localized regime, whereas they show ordi- 
nary localization for short-range disorder that mixes the 
valleys and restores TRS In contrast, armchair GNRs 
were generally considered to be i n the orthogonal sym- 
metry class corresponding to TRSP32H 

In this paper, we investigate carefully the symmetries 
of graphene nanoribbons and their effect on the conduc- 
tance in the strongly localized and diffusive regime when 
only long-range disorder (that does not mix the valleys) 
is present. In particular, we show that in contrast to 
common belief the symmetry classification of armchair 
GNRs depends on whether they are metallic or semi- 
conducting. While semiconducting armchair GNRs are 
found to be in the orthogonal symmetry class, metallic 
armchair GNRs exhibit a hidden pseudovalley structure 
which leads to effective TRS-breaking and places metal- 
lic armchair GNRs into the unitary symmetry class. This 
pseudovalley structure also leads to a perfectly conduct- 
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ing channel in metallic armchair GNRs. 

In addition, we show that zigzag GNRs can exhibit 
an unexpected and strong source of intervalley scatter- 
ing, even for long-range, smooth potentials. When the 
magnitude of the disorder potential exceeds the Fermi 
energy, electron-hole puddles are formed and valley scat- 
tering can be mediated by the edge state. In this case, 
TRS is restored in zigzag GNRs and the PCC vanishes. 
This puts an additional restriction on the disorder po- 
tential (apart from being long-ranged) in order to realize 
the unitary symmetry class in zigzag GNRs. 

Weak localization effects in the conductance in 
the diffusive regime of extended bulk graphene 
have been st udied extensively both theoretically and 
experimentallyPil2iH!U i n contrast, we are aware only 
of a single theoretical work for (quasi-onedimensional) 
GNRs that finds weak localization behaviorP^l Here we 
study the quantum transport properties of GNRs in the 
diffusive regime systematically, and show that their be- 
havior is in agreement with their symmetry classifica- 
tions. In particular, we find suppression of weak local- 
ization in zigzag GNRs without intervalley scattering and 
metallic armchair GNRs. In contrast, zigzag GNRs with 
strong intervalley scattering and semiconducting arm- 
chair GNRs exhibit weak localization behavior. The sym- 
metry classification also reflects itself in the conductance 
fluctuations of the GNRs. 

The paper is organized as follows: In Sec. [II] we intro- 
duce the tight-binding and Dirac Hamiltonian describ- 
ing the electronic structure of graphene, and briefly de- 
scribe our numerical method. We investigate the valley 
scattering properties of a long-range disorder potential in 



Sec. Ill showing that there is an unexpected source of 
valley scattering for zigzag GNRs. In Sec. |IV| we classify 
the different types of GNRs according to their symme- 
try and study their quantum transport properties. We 
finally conclude in Sec. [Vj 



II. GRAPHENE HAMILTONIAN AND 
QUANTUM TRANSPORT 

A. Hamiltonian 





FIG. 1: (a) The graphene honeycomb lattice. The A and B 
sublattices are indicated as solid and open dots, respectively, 
(b) zigzag and (c) armchair graphene nanoribbons. The width 
W of GNRs is measured between the first rows of missing 
atoms (shown in grey). 

For a sufficiently smooth potential V(x) and in the 
low-energy limit, excitations with energy e obey the Dirac 
equation 



where the Hamiltonian 



H = v f t ® {a ■ p) + V(x) t ® cr 
acts on a four-component spinor wave function 



(2) 



(3) 



(4) 



The Hamiltonian is written in the valley-isotropic form 
introduced in Ref. 1311 t% and <7j denote the Pauli matrices 
in valley and sublattice space, respectively (to and cro are 
the respective unit matrices), and \& j and ^ with j € 
{^4, B} are the wave function amplitudes on the different 
sublattices in the K and X'-valley. The Fermi velocity is 
denoted as v-p and p = —ih(d xi d y ) is the two-dimensional 
momentum operator, with the orientation of x and y-axis 
as indicated in Fig. [lja). A magnetic field is included 
through the minimal coupling p — > p + eA(x) with — e 
the electron charge. 



We describe the electronic structure of graphene using 
a tight-binding model for the honeycomb lattice [shown 
inFig.Qa)], 

with one orbital \i) per atom and constant hopping tij = t 
only between nearest neighbors. We allow for an on-site 
potential V{x) that is evaluated at the positions of the 
carbon atoms Xi. A magnetic field is included through 
the substitution £y — > t x exp[if J * dxA(x)], where 

A(x) is the magnetic vector potential. 



B. Numerical quantum transport in the 
tight-binding model 

To support our analytical predictions, we perform nu- 
merical computations of the quantum transport proper- 
ties of graphene nanoribbons cut out of the graphene lat- 
tice [examples of zigzag and armchair graphene nanorib- 
bons are shown in Fig.HTb) and (c)]. 

A potential term V(x) is only introduced in a finite 
part of the system (the scattering region) , the remaining 
parts, i.e. perfect semi-infinite nanoribbons then serve 
as leads (with a Fermi energy identical to the scatter- 
ing region). We compute the lattice Green's function of 
the system using an adaptive recursive Green's function 
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FIG. 2: (Color online) (a) Example of a realization of the 
impurity potential Q. (b) Numerically computed probability 
distribution of the value of the potential V(x) on a given point 
x. For both (a) and (b), pi mp = 0.05, £ = 2a, and 5 = 0.05 1. 



technique.^ From the Green's function we then obtain 
the scattering matrix using the Fisher-Lee relation for 
tight-binding systems.^ In particular, we compute the 
amplitudes t n>m for transmission from mode m to mode n 
between two leads, and the amplitudes r„ jI7l for reflection 
from mode m to mode n in the same lead. The electri- 
cal conductance G is th en obtained using the Landauer- 
Biittiker formalisrrPEH, 



G 



(5) 



where Go = 2t- is the conductance quantum including 
the spin degree of freedom. 



III. INTER- VALLEY SCATTERING IN 
DISORDERED GRAPHENE 

A. Impurity potential 

We use a model for an impurity potential that is com- 
monly used in the study of disordered graphene (e.g. 
Refs. 27 36 39). The potential is assumed to consist of a 
set of independent impurities with a Gaussian potential 
profile: 



A'i 



V(x) = Sj exp 



2e 



(6) 



Here Xj is the position of the j'-th scattering center. We 
use Ni = p- lmp iV a randomly distributed Gaussian scat- 
terers, where N a denotes the total number of lattice sites 
in the disordered region, and Pi mp < 1 is a constant that 
determines the relative amount of scatterers. Further, we 
choose the impurity strength 8j randomly from the inter- 
val [—5, 5] and use a constant range £ for all impurities. 
Figure [2]ja) shows an example of the potential landscape 
of this type of impurity potential. 

The potential ^ can describe both short-range and 
long-range impurities by varying the impurity range 
£. For £ > a, the resulting potential varies smoothly 



on the lattice scale, and the Dirac Hamiltonian ([3| is 
applicable.^ 

The magnitude of the fluctuating impurity potential is 
best described by its root mean square value. For £ 3> a 
(in practice it is enough to have £ > a) it is given as 



I 4-7T Pi 



3^3 



imp S 1 
a 



(7) 



where (...) denotes an impurity average. Figure [2] shows 
the probability distribution P (V(x)) for finding a partic- 
ular potential value V{x) at a point x. The distribution 
is Gaussian-like, with a width given by ay. 

In the limit of Af 3> £ > a (with Af the Fermi wave 
length), the tra nsport mean free path is given in Born 
approximation a d 26 * 36 * 
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with the dimcnsionless correlator 
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B. Bulk graphene versus zigzag nanoribbons 

The parameter £ of the impurity potential ^ deter- 
mines the smoothness of the potential. It is generally 
accepted that for £ > a there is only little intervalley 
scattering. However, the evidence for this was always 
only indirect [ 23 l 27 l 36 l 37 ] an d no quantitative discussion of 
intervalley scattering does exist for this type of poten- 
tial. Recently, it has only been attempted to quantify 
the intervalley scattering for short-range lattice defects™ 
Since the knowledge of the degree of intervalley scattering 
will be important in the following section, we first investi- 
gate the intervalley scattering for the impurity potential 
Our findings show that caution must be exerted, as 
the presence of zigzag edges may lead to enhanced in- 
tervalley scattering even if the impurity potential is very 
smooth. 

We can numerically measure the intervalley scattering 
if we consider a wire along the x-direction, using either 
periodic boundary conditions in y-direction (making the 
system equivalent to an armchair carbon nanotube) or 
zigzag boundaries. In both cases the valleys K and K' 
project onto different values of the longitudinal momen- 
tum k x . The scattering states in the leads have a definite 
Bloch momentum k x , and a mode m can thus be uniquely 
assigned to a valley (this is not possible for armchair rib- 
bons, where the two valleys project on the same momen- 
tum). It is then possible to decompose the numerically 
computed transmission and reflection probabilities into 
an intravalley and intervalley part. The total probability 
of reflection into the other valley is given as 



R 



inter 



E 



(9) 
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FIG. 3: (Color online) Probability for intervalley scattering 
in reflection [Eq. |9|] as a function of impurity strength 5 
for a ribbon with periodic boundary conditions (left panel) 
and zigzag boundaries (right panel). The ribbons have width 
W ~ 500a and length L ~ 2500a, and data is shown for 
E = 0.017* (black solid lines) and E = O.lt (red solid lines). 
For comparison, the right panel (with zigzag boundaries) also 
contains the data for E = O.lt with periodic boundary condi- 
tions as a red dashed line. The data was obtained by averag- 
ing over 100 different impurity configurations with potential 
parameters pi mp = 0.05 and £ = 2a. 



junctions (in the bulk, but also at the zigzag edge), when 
its magnitude becomes larger than the Fermi energy E: 
For E = O.lt and small 8, the intervalley scattering for 
the zigzag case follows the result obtained with periodic 
boundary conditions (shown for comparison as a dashed 
line in the right panel of Fig. [3]). For larger 8, when the 
maxima of the potential become comparable or greater 
to E = O.lt (which is already the case around 8 = 0.05 1 
as seen from Fig. [2]), the intervalley scattering rate raises 
rapidly and reaches up to 10%. For smaller Fermi en- 
ergy E = 0.017i, this regime is reached already for much 
smaller 8, and we find perfect valley mixing close to 50%. 

Hence, although the impurity potential itself indeed 
does not scatter valleys significantly if £ > a, caution 
must be exerted if zigzag edges are present. Intervalley 
scattering can be very large if the impurity potential lo- 
cally crosses the Dirac point and hence locally forms pn- 
junctions. We believe that this unexpected valley scat- 
tering may also explain many not understood numerical 
findings of the past. 



IV. SYMMETRY AND THE CONDUCTANCE 
OF GRAPHENE NANORIBBONS 

A. Symmetries and quantum transport in 
disordered wires 



while the total reflection probability is given as R — 
J2nm \ r n,m\ 2 - The intervalley transmission probability 
can be defined analogously. 

In Fig. [3] we show the fraction of intervalley scatter- 
ing in the reflection probability, Ri n t OT /R, as a function 
of the impurity strength 8. The rationale for measuring 
the intervalley scattering in the reflection is that trans- 
mission may contain a sizeable part of ballistic processes 
without any scattering; in contrast, reflection only occurs 
after at least one scattering event. For periodic bound- 
ary conditions and £ = 2a we indeed find only very little 
intervalley scattering (left panel of Fig. |3j. The inter- 
valley scattering rises with increasing impurity strength 
8. This is to be expected, as for fixed £ the potential 
becomes steeper as 8 is increased, and hence interval- 
ley scattering becomes more likely. Nevertheless, for the 
given parameters, the fraction of intervalley scattering 
remained below 10~ 3 . 

However, we obtain a very different picture for zigzag 
graphene nanoribbons (right panel of Fig. [3]). Although 
we use the same impurity potential as in the case of peri- 
odic boundary conditions, we find an intervalley scatter- 
ing that is 3 orders of magnitude larger. The only obvious 
difference is the presence of the zigzag edge. Indeed, it 
was shown that even a smooth pn-j unction in a zigzag 
nanoribbon strongly scatters valleys, as both valleys are 
connected by the edge stateP^ We believe that our nu- 
merical findings can be explained fully by the fact that 
the impurity potential locally leads to many smooth pn- 



The quantum transport properties of disordered quan- 
tum wires are universal and determined by their symme- 
tries only.Esl I n particular, the symmetry class is deter- 
mined by the presence or absence of time-reversal sym- 
metry (TRS) T, which is an antiunitary symmetry. A 
system may belong to one of the three Wigner symme- 
try classes: unitary if TRS is broken, and orthogonal or 
symplectic if TRS is present with T 2 = +1 or T 2 = — 1, 
respectively. A system with TRS obeys THT^ 1 = H . 

Note that a Hamiltonian may possess a TRS that is, 
however, irrelevant: This is the case if the Hamiltonian 
decomposes into independent blocks and the symmetry 
connects only between them. In this case the TRS has 
no influence on the quantum transport properties (except 
guaranteeing a degeneracy between the blocks), and the 
symmetry class is determined by intrablock symmetries 
only. Below, we identify the appropriate symmetries for 
the case of graphene nanoribbons. 



B. Symmetries of graphene nanoribbons within the 
Dirac approximation 

1. Bulk symmetries of graphene and boundary conditions 

The bulk Dirac Hamiltonian ^ commutes with four 
antiunitary symmetrieflSMMEIl 

% — Ti®a y C for i e 0, x, y, z, (10) 
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where C denotes complex conjugation. Each of these an- 
tiunitary symmetries can play the role of an (effective) 
TRS. 

It is easy to see that T y 2 = 1, whereas 7~ 2 = — 1 for 
i £ {0,x,z}. T y represents the (true) TRS that connects 
the two valleys. T x is the valley symmetry that guar- 
antees the Kramer's degeneracy of both valleys (since 
T x — — 1). It should be noted that 7o and T z only differ 
by a phase in the two valleys and are thus equivalent. 
In fact, one can write down a whole family of equivalent 
effective intravalley TRSfP 



7oz($) = (cos To + ism-dT z ) 



OyC . 



(11) 



From this family, a single antiunitary symmetry will sur- 
vive in the case of metallic armchair nanoribbons, as we 
will show below. The presence of a magnetic field breaks 
all four symmetries 71 . 

The other important symmetries of graphenc, chiral 
symmetry (r z ® <j z )H(t z ®a z ) = —H, and particle-hole 
symmetry (t <8> <j x )H*(t ® a x ) = —H, are broken by 
the potential term V(x). Hence, we do not expect to see 
universality classes beyond the three Wigner classes.^ 

It is well-known that bulk graphene with long-r ange 
scatterers belongs to the symplectic symmetry class 122125] 
In this case valleys are not mixed and the true TRS T y 
as well as the valley symmetry T x are irrelevant, and the 
effective intravalley TRS 7o, z determines the symmetry 
class. 

In a graphene nanoribbon, the antiunitary symmetries 
71 must also be compatible with the boundary conditions. 
Boundary conditions in th e Dirac equation can be written 
generally in the forrrPUl 

^>{x) = Mb^(x) for x on boundary 6, (12) 

where M. b is a Hermitian 4 x 4- matrix. A graphene 
nanoribbon is then symmetric with respect to %, if 



[M b ,Ti] = 0, 



(13) 



where [A, B) = AB—BA denotes the commutator. Below 
we now specialize to the cases of zigzag and armchair 
nanoribbons. 



2. Zigzag nanoribbons 



The boundary condition for a zigzag GNR reai 



±tv ® o 7 for b = 1,2. 



(14) 



This boundary condition does not mix valleys which thus 
remain a good quantum number. The system is symmet- 
ric with respect to T y and 7~ x , but the boundary condi- 
tion breaks 7o,z- However, since valleys are not mixed for 
long-range impurities, the valley-offdiagonal symmetries 
y are not relevant. Since all intravalley TRS are broken 
by the boundary conditions, a zigzag graphene nanorib- 
bon resides in the unitary symmetry class, as shown in 
Ref . m\ 



These considerations hold as long as there is no inter- 
valley scattering due to local pn-junctions at the zigzag 
edge. If there is, the valleys are strongly mixed, and the 
only remaining symmetry is the TRS of the tight-binding 
lattice (7tb = C, 7^, = 1). The Dirac equation cannot 
capture scattering between the valleys via the edge state 
since a continuous equation cannot represent the finite 
size of the zigzag GNR Brillouin zone that is at the heart 
of the scattering mechanism.^ The zigzag GNR is then 
in the orthogonal symmetry class. 



3. Metallic armchair nanoribbons 



The boundary condition for an armchair GNR reads^Sl 



Mr, 



t ■ <J 



v>b = (sin cos$b, 0) (15) 



with b = 1,2, and t points in the direction of the GNR. 
This boundary condition strongly mixes valleys, and the 
relative valley-angle $2 — $1 between the two boundaries 
of the armchair GNR depends on the width of the rib- 
bon. For example, for an armchair GNR in ^-direction 
as shown in Fig. [ljc), $1, = —2Kx\, where K = 47r/3aP^ 
A metallic armchair nanoribbon has $1 = #2 = i-e. 



Mx=M 2 =M. 



(16) 



In this case, u ■ t commutes with both the Hamiltonian 
(including disorder) and the boundary condition, and we 
may choose the solutions of the Dirac equation as eigen- 
states of ut. The solutions can thus be grouped into two 
new pseudov alleys, K-jz and Kt^, that remain uncoupled 
for long-range disorder. The pseudovalley description is 
obtained from the usual valleys by means of the rotation 



n 



-i-KT-e/A — 



(17) 



The metallic armchair boundary condition ( 16 ) is sym 



metric with respect to T y and 7o z (^) with d equal to 
the valley-angle of the boundary condition, whereas T x 
is broken. In the pseudovalley space they take the form 



= r v 
7») = 



CTyC , 



OyC . 



(18) 



where T 7 ^ = IZTTZ^ ■ Hence, both symmetries are com- 
pletely offdiagonal in valley space and are not relevant 
for determining the symmetry class in the case of long- 
range potential. In the absence of any intra-pseudovalley 
TRS, metallic armchair GNRs also belong to the unitary 
symmetry class. Note that Ref. already discussed ef- 
fective TRS breaking in the context of the lowest mode of 
metallic armchair GNRs in a ring geometry. Our analy- 
sis shows the more general result that, in the low-energy 
limit, metallic armchair GNRs belong to the unitary sym- 
metry class without a special geometry and regardless of 
the number of modes. We will confirm this using numer- 
ical simulations below. 



6 




-2 2 

kW/n 



-2 2 
kW/n 



-2 2 

Win 



FIG. 4: (Color online) Pseudovalley resolved band struc- 
ture of a metallic armchair nanoribbon after rotation 1Z. see 



Eq. (171 



A metallic armchair GNR features a band structure 
e(fc) (with fc the Bloch wave vector) with two non- 
degenerate gapless linearly dispersing bands and pairs of 
two- fold degenerate hyperbolic bands (shown in Fig. Ej) P 
We can use the conserved antiunitary symmetries to un- 
ravel the pseudovalley structure: Both Tj 1 and 7 u J(i9) 
lead tcP 



(19) 



Hence, the two counter-propagating gapless linear modes 
are (Kramer's) partners residing in different pseudoval- 
leys. In addition, every pseudovalley contains a set of 
non-degenerate hyperbolic bands, as shown in Fig. |4j 

The pseudovalley structure also reveals that the metal- 
licity (i.e. absence of a gap) of metallic armchair nanorib- 
bons is of topological origin: In pseudovalley space, the 
boundary condition reads 



M 



t ■ <j . 



(20) 



and thus take s the form of infinite mass boundary 
conditions, 6 19 with an infinite mass of opposite sign on 
the opposite edges (if both edges had the same sign of 
mass, A4i = — M 2 )- Hence, the metallic armchair GNR 
effectively exhibits a domain wall with a sign change 
in m ass an d thus supports a gapless linearly dispersing 
mode™ 

Finally, it must be emphasized that the pseudovalley 
structure is only valid for energies close to the Dirac 
point. For higher energies trigonal w arpin g breaks the 
symmetry between the K and K' valleyPSESlt introduces 
terms proportional to r z in the Hamiltonian which then 
does not commute any more with v ■ t and thus invali- 
dates the notion of pseudovalleys. In fact, the degeneracy 
of hyperbolic bands (a consequence of the pseudovalley 
structure) in the tight-binding model is only true close to 
the Dirac point, as shown in Fig. [5]ja). For large energies 
it is absent [Fig. (5jc)]. 



Jf.. Semiconducting armchair nanoribbons 



An armchair GNR is semiconducting, if $1 ^ $2 in 
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FIG. 5: (a), (b) Tight-binding band structures of wide arm- 
chair graphene nanoribbons close to the Dirac point, (a) 
Metallic armchair GNR with W = 501 a. The non-linear 
bands are approximately degenerate for low energies, (b) 
Semiconducting armchair GNR with W = 502 a. (c) Full 
tight-binding band structure of a narrow metallic armchair 
GNR with W — 7.5 a. Clearly the band degeneracy is lifted, 
(d) Narrow zigzag GNR with W = 19/V3a. 



ary conditions, 



Mi + M 2 - 



(21) 



In this case, it is not possible to separate any valley struc- 
ture, and hence there is no degeneracy of bands even close 
to the Dirac point [Fig. [5][b)]. T y is the only symme- 
try of the problem, and hence semiconducting armchair 
nanoribbons with long-range disorder belong to the or- 
thogonal symmetry class. 

Our findings for the different types of GNRs are sum- 
marized in Table lU 



C. Perfectly conducting channels 

1. Previous work 

One of the most striking features of a zigzag GNR in 
the absence of intervalley scattering is the presence of a 
perfectly conducting channel (PCC)P^In this case, one 
of the transmission eigenvalues is exactly one, such that 
G/Gq > 1, regardless of the strength of the disorder. 
As explained in Ref. 1231 within a single valley, a zigzag 
GNR has unequal numbers p, q of right-moving and left- 
moving channels, respectively. This can only occur in the 
absence of TRS and limits the conductance from below 



Eq. (151, i.e. if the two boundaries have different bound- 



G/G >\ P -q\ 



(22) 
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GNR 


symmetry class 


metallic armchair 


unitary 


zigzag, no intervalley scattering due 
to pn-j unctions at edge 


unitary 


zigzag, intervalley scattering due to 
pn-junctions at edge 


orthogonal 


semiconducting armchair 


orthogonal 



a) 3.0 



W= 251a (metallic) 
W— 252 a (semicond.) 



TABLE I: Summary of symmetry classification of GNRs with 
long-range disorder and energies close to the Dirac point. For 
zigzag GNRs we must distinguish whether there is intervalley 
scattering due to local pn-junctions at the zigzag edge. 



In a zigzag GNR the difference in right- and left-movers, 
\p — q\ = 1, is due to the zigzag edge state that connects 
the two valleys [see Fig. [5^d)]. 



2. Armchair nanoribbons 

From our previous analysis, we found that metallic 
armchair GNRs also exhibit broken (effective) TRS and 
have \p — q\ = 1 within a single pseudovalley (Fig. [4]). 
The inequality of left and right-movers in this situation 
is associated with the linearly dispersing modes of op- 
posite velocity in the two pseudovalleys. Hence, we also 
expect a PCC in this situation. In contrast, a semicon- 
ducting armchair GNR belongs to the orthogonal group 
and we expect to see ordinary localization. 

This is confirmed by our numerical simulations in 
Fig. [6^ a) : We clearly see different localization behavior 
for semiconducting and metallic armchair GNRs, with 
the latter saturating at G/Gq = 1 for long ribbons, and 
thus exhibiting a PCC. It should be emphasized that 
this PCC is not identical to a single channel of the GNR 
(i.e. the edge state in a zigzag GNR or the linearly dis- 
persing mode in a metallic armchair GNR), instead the 
unit transmission eigenvalue corresponds to a superposi- 
tion of all channels. 

A previous numerical study^ of metallic armchair 
GNRs with long-range disorder found an only approxi- 
mately unit conductance quantization in the single-mode 
regime, and the absence of a PCC in the multi-channel 
regime. In contrast, our numerical data shows a PCC 
even in the multi-mode regime [three open channels in 
Fig. [6ja,)]. This apparent contradiction is resolved by 
noting that the simulations of Ref. [55] used very narrow 
GNRs (W — 7.5a s» 1.8 nm), where the second channel 
only opens at energies far from the Dirac point [around 
0.3i in Fig. gc)]. For those high energies, trigonal warp- 
ing breaks the pseudovalley structure and the PCC van- 
ishes. 



3. Zigzag nanoribbons 

In our simulations we have carefully chosen the param- 
eters of the disorder potential such that bulk valley scat- 




3 4 5 
Lx 10 3 a 



FIG. 6: (Color online) Average conductance of (a) armchair 
and (b) zigzag GNRs as a function of ribbon length. In all 
systems the Fermi energy corresponds to three open channels, 
(a) Semiconducting armchair GNR (black) with W = 251 a 
and metallic armchair GNR (blue) with W = 252 a. In both 
cases, Fermi energy E — 0.0174 and the disorder parameters 
are p imp = 0.05, £ = 2.0 a and 8 = 0.08 4. (b) Zigzag GNR 
with different widths and disorder parameters (with pi mp = 
0.05 and £ = 2.0a fixed). Black: W = 434.5/^0 w 251a, 
E = 0.022t, 8 = 0.08 1. Red: W = 88/V3a « 51a, E = 
0.11 1, 8 = 0.084. Orange: W = 88/V3a w 51a, E = 0.11 t, 
8 = 0.04 1. Green (inset): W = U.5/V3a « 26a, E = 0.22t, 
8 — 0.08 1. The data in (a) and (b) was averaged over 50-200 
impurity configurations. 



tering is indeed negligible, as shown in Sec. |IIIB| How- 
ever, there we also observed that valley scattering could 
be large in the presence of zigzag edges. This has also 
consequences for the PCC in zigzag GNRs: In Fig. |6^b) 
we show the average conductance as a function of rib- 
bon length for various energies, disorder parameters and 
width. In particular, zigzag GNRs do not show a PCC 
for a disorder where a metallic armchair GNR very well 
did (black line in Fig. l^b) and blue line in Fig. |6][a), 
respectively). Only if the amplitude of the disorder po- 
tential is smaller than the Fermi energy, i.e. if no local 
p-n junctions are formed, a PCC can be observed [orange 
line in Fig.^b)]. However, as the disorder has to be cho- 
sen weaker, the conductance saturates only for very long 
GNRs. This breakdown of the PCC due to valley scat- 
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tering mediated through the zigzag edge state has not 
been observed in previous studies that dealt with narrow 
ribbons (W ~ 5 a « 1.2 nm in Ref. 23) at energies further 
away from the Dirac point [for an example of a PCC in 
this case, see inset of Fig. |6jb)]. 

4- Zigzag vs. metallic armchair nanoribbons 

In summary, zigzag and metallic armchair GNRs both 
exhibit a PCC. In metallic armchair GNRs, its observa- 
tions requires small energies close to the Dirac point, it 
vanishes when the Fermi energy is larger and in a regime 
where trigonal warping becomes effective. In contrast, 
for zigzag GNRs the Fermi energy must be larger than 
the potential amplitude; otherwise the valleys are cou- 
pled and the PCC vanishes. Hence, in metallic armchair 
GNRs the Fermi energy should be smaller than an in- 
trinsic energy scale (trigonal warping) , whereas in zigzag 
GNRs the Fermi energy should be larger than an extrin- 
sic energy scale (disorder potential) in order to observe 
the effective TRS-breaking and a PCC. 



D. Magnetoconductance 




B/T 

FIG. 7: Normalized magnetoconductance of disordered 
GNRs: semiconducting armchair GNR (black) with W = 
502a, and E = 0.017i (six open channels), metallic armchair 
GNR (blue) with W = 501a, and E = 0.017* (7 open chan- 
nels), and zigzag GNR with W = 865.5/v^a » 500 a and 
E — 0.017 1 (red, 5 open channels) and E = 0.1 1 (violet, 35 
open channels). In all cases, the GNR length was L « 2500 a, 
and the disorder parameters Pi mp = 0.05, £ = 2.0a and 
5 — 0.08 1 [as in Fig. [6ja)] . The data was averaged over 600 
impurity configurations. 



Apart from the PCC that reveals itself mainly in the 
strongly localized regime, the symmetry class also influ- 
ences the conductance in the diffusive regime. Due to 
quantum-coherence corrections, the conductance can be 
either smaller (weak localization, WL), larger (weak an- 
tilocalization, WAL) or equal to the classical conductance 
for the orthogonal, symplectic, and unitary symmetry 
classes, respectively. 18 These quantum coherence correc- 
tions reveal themselves in the magnetoconductance, in 
particular in the change of the disorder averaged con- 
ductance with magnetic field: 

(5G(B)) = {G(B))-{G(B = 0)), (23) 

where B is a magnetic field perpendicular to the GNR 
and G(B) the conductance for given field B. Since a 
magnetic field breaks TRS, every GNR is in the uni- 
tary symmetry class for large enough B. Hence, for 
large enough B, (SG(B)) > (WL) in the orthogonal 
class, (SG(B)) < (WAL) in the symplectic class, and 
(SG(B)) = (suppressed WL) in the unitary class. 

Figure [JJ summarizes our results of magnetotransport 
simulations for zigzag and armchair GNRs. As expected 
from our symmetry considerations in Sec. |IVB| we ob- 
serve WL behavior for semiconducting armchair GNRs 
(black curve), whereas the WL correction is suppressed 
in the metallic armchair GNR (blue curve). The results 
for the zigzag GNRs demonstrate again the importance of 
intervalley scattering: For small Fermi energy (red line) 
the intervalley scattering is large [parameters as for the 
black line in Fig.[3^b)] leaving only the tight-binding TRS 
7tb- Thus we observe WL behavior just as in the case of 



the semiconducting armchair GNR, as both belong to 
the orthogonal symmetry class. Only if valley scattering 
is suppressed for larger energies (violet curve), we also 
observe a suppression of the WL correction, as expected 
from the unitary symmetry class (there is some residual 
intervalley scattering in this case, preventing complete 
suppression as in the m etallic armchair GNR). 

Previous studies^ 21221 assum ed that the role of edges is 
only to introduce valley scattering and hence expected 
WL behavior in GNRs. In contrast, our study has shown 
that the type of edge, and even the distance between 
opposite edges is crucial to understand the magnetocon- 
ductance of GNRs in the quantum regime. 

Random matrix theory^ (RMT) predicts a universal 
value of the quantum-coherence correction in the limit of 
a large channel number. This value only depends on how 
far the system has approached the diffusive limit mea- 
sured by the parameter s — 2L/7rZ tr . From Eq. pi we 
obtain s w 1.6 for the chosen disorder parameters [see 
caption of Fig. [JJ , which agrees well with the value es- 
timated from the average conductance in the simulation 
(s « 2). The value of the WL correction in the orthogonal 
symmetry class is therP 1 lim B ^ 00 ((5G'(S)/Go) ~ 0.26. 
The WL correction obtained from our numerical simula- 
tions agrees reasonably with this prediction, given that 
the number of channels is still small. 



E. Universal conductance fluctuations 

In addition to the quantum-coherence correction to 
the (average) conductance, RMT also predicts universal 
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— zigzag, E = O.lf 

— semiconducting ac 

— zigzag, E = 0.017r 
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FIG. 8: Universal conductance fluctuations as a function 
of magnetic field for semiconducting armchair (ac) GNRs 
(black), metallic armchair GNRs (blue), and zigzag GNRs 
(red: E — 0.017i, large intervalley scattering; violet: E = 
O.lt, small intervalley scattering). The dashed lines show the 
RMT predictions for the value of the UCFs for s = 1.6. Pa- 
rameters as in Fig. [7] 



values of the conductance fluctuations Var(G/Go) (uni- 
versal conductance fluctuations, UCF).^ Both, zigzag 
and metallic armchair GNRs are in the unitary sym- 
metry class but have degenerate (pseudo) valleys. As a 
consequence, we expect UCFs of four times the univer- 
sal value of the unitary class. In magnetic field, the 
(pseudo) valleys remain independent, but their degener- 
acy is broken, since all of the antiunitary symmetries 
( 10 1 are broken by magnetic field. Hence the UCFs take 
twice the value of the unitary class. A semiconducting 



armchair GNR does not allow for a decomposition into 
independent blocks and thus the UCFs take the value of 
the orthogonal and unitary class in the absence and pres- 
ence of a magnetic field, respectively. Since the UCFs in 
the orthogonal symmetry class are twice as large as in the 
unitary symmetry class, the UCFs in zigzag and metallic 
armchair GNRs are twice as large as for semiconducting 
armchair nanoribbons, both in the absence and presence 
of a magnetic field. 

Fig. [8] shows the universal conductance fluctuations 
as obtained from our numerical simulations. We in- 
deed observe that the UCFs of metallic armchair GNRs 
and zigzag GNRs with little intervalley scattering are al- 
ways significantly larger than the UCFs of semiconduct- 
ing armchair GNRs and zigzag GNRs with large inter- 
valley scattering. We even find good quantitative agree- 
ment with the RMT values. Only the UCFs for semi- 
conducting armchair GNRs and zigzag GNRs with large 
intervalley scattering are somewhat larger than the theo- 
retical prediction, probably due to the still small number 
of channels. 



V. CONCLUSIONS 

We have carefully investigated the symmetry classifica- 
tions of graphene nanoribbons with long-range disorder 
in the Dirac limit, and studied their imprints on the con- 
ductance. Table |T] summarizes our findings. 

In principle, all of the considered graphene nanorib- 
bons are time-reversal symmetric in the absence of a mag- 
netic field. However, if intervalley scattering is absent 
(hence the condition of long-range disorder), this true 
TRS that connects the two valleys is irrelevant, and the 
type of boundary is decisive for the symmetry properties. 

In particular, we have found that in the case of arm- 
chair GNRs (that up to now were generally assumed to 
be in the orthogonal symmetry class) it is necessary to 
distinguish between semiconducting and metallic vari- 
ants: While semiconducting armchair GNRs inevitably 
mix valleys and thus belong to the orthogonal symmetry 
class, metallic armchair GNRs have a hidden pseudoval- 
ley structure that together with the boundary conditions 
places them into the unitary symmetry class. 

Zigzag graphene nanoribbons have already 
previously^! been identified to belong to the uni- 
tary symmetry class. However, we have shown for this 
classification it is necessary that the Fermi energy is 
larger than the disorder potential fluctuations. Oth- 
erwise, local pn-j unctions at the zigzag edge act as 
strong intervalley scatterersP^ We have demonstrated 
numerically that the intervalley scattering due to this 
mechanism can lead to complete valley mixing, although 
the disorder potential alone would not scatter valleys. 
Hence, for zigzag nanoribbons to be in the unitary class 
it is not enough to be in the Dirac limit and to have 
long-range disorder, there is also a restriction on the 
magnitude of the potential with respect to the Fermi 
energy. 

The symmetries of the GNRs also have a strong influ- 
ence on their quantum transport properties. In a metal- 
lic armchair GNR, the pseudovalley structure manifests 
itself most conspicuously in a perfectly conducting chan- 
nel, i.e. a lower bound of one conductance quantum even 
in a strongly disordered nanoribbon. 

The perfectly conducting channel reveals itself most 
clearly in the strongly localized regime, but the symme- 
tries of the GNRs also manifest themselves in the diffu- 
sive regime. We showed that weak localization is strongly 
suppressed in metallic armchair GNRs as well as zigzag 
GNRs with little intervalley scattering. In contrast, semi- 
conducting armchair GNRs and zigzag GNRs with a dis- 
order potential amplitude larger than the Fermi energy 
exhibit weak localization (instead of weak antilocaliza- 
tion expected for bulk graphene with long-range disor- 
der) due to intervalley scattering at the armchair edges 
and local pn-j unction at the zigzag edge, respectively. 

In addition, the interplay of symmetry classes and the 
degeneracy of the (pseudo)valley structure of metallic 
armchair GNRs and zigzag GNRs with little intervalley 
scattering leads to larger conductance fluctuations than 
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in semiconducting armchair GNRs and zigzag GNRs with 
intervalley scattering. 
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